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' Asymptotic expansions of Green functions and spectral densities associ- 

ated with partial differential operators are widely applied in quantum field 
theory and elsewhere. The mathematical properties of these expansions can 
be clarified and more precisely determined by means of tools from distribu- 
Q | tion theory and summability theory. (These are the same, insofar as recently 

the classic Cesaro-Riesz theory of summability of series and integrals has 
been given a distributional interpretation.) When applied to the spectral 
analysis of Green functions (which are then to be expanded as series in a pa- 
rameter, usually the time), these methods show: (1) The "local" or "global" 
dependence of the expansion coefficients on the background geometry, etc., is 
determined by the regularity of the asymptotic expansion of the integrand at 
the origin (in "frequency space"); this marks the difference between a heat 
kernel and a Wightman two-point function, for instance. (2) The behavior 
of the integrand at infinity determines whether the expansion of the Green 
function is genuinely asymptotic in the literal, pointwise sense, or is merely 
valid in a distributional (Cesaro- averaged) sense; this is the difference between 
the heat kernel and the Schrodinger kernel. (3) The high-frequency expan- 
sion of the spectral density itself is local in a distributional sense (but not 
pointwise). These observations make rigorous sense out of calculations in the 
physics literature that are sometimes dismissed as merely formal. 
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I. INTRODUCTION 



The aim of this article is to study several issues related to the small-t behavior of var- 
ious Green functions G(t,x,y) associated to an elliptic differential operator H. These are 
the integral kernels of operator-valued functions of H, such as the heat operator e~ tH , the 
Schrodinger propagator e~ ltH , various wave-equation operators such as cos(ty/~H), the oper- 
ator e~ is ^ that solves a certain elliptic boundary-value problem involving H, etc. All these 
kernels are expressed (possibly after some redefinitions of variables) in the form 

/■oo 

G(t,x,y) = g(t\)dE x (x,y), (1) 
Jo 

where E\ is the spectral decomposition of H, and g is a smooth function on (0, oo). 

Each such Green function raises a set of interrelated questions, which are illumined by 
the simple examples listed in the Appendix: 

(i) Does G(t.x,y) have an asymptotic expansion as t 10? For the heat problem, ( |A1|) , 
it is well known0i that 

oo 

K(t, x, x) ~ (Anty d/2 a n (x, x)t n/2 , (2a) 

n=0 

where d is the dimension of the manifold M. and ao(x,x) = 1. Similar formulas hold off- 
diagonal; for example, if Ai C R d and the leading term in H is the Laplacian, then 

oo 

K(t, x, y) ~ (47rt)- d / 2 e-l^l 2 / 4t £ a n (x, y)t n ' 2 . (2b) 

n=0 

In the case ( |A7b| ), the elementary heat kernel on R , all a n = except the first. In fact, 
this is true also of ( |Allb| ), the elementary Dirichlet heat kernel on (0, 7r), because as t goes 
to the ratio of any other term to the largest term (e~^ x ~ y ^ / 4t ) vanishes faster than any 
power of t. In particular, therefore, the expansion (fj) for fixed (x, y) G (0, n) x (0, it) does 
not distinguish between the finite region (0, 7r) and the infinite region R. (However, the 
smallness of the two nearest image terms in ( |Allb| ) is not uniform near the boundary, and 
hence f£ K(t,x,x) has an asymptotic expansion (47rt) -1 / 2 Yl™=o A n with nontrivial higher- 
order terms A n .) This "locality" property will concern us again in questions (iv) and (v). 

The Schrodinger problem, (|A2|) , gives rise to an expansion @ that is formally identical 
to (0) (more precisely, obtained from it by the obvious analytic continuation) .i'i However, 
it is obvious from ( |A12b| ) that this expansion (which again reduces to a single term in the 



examples ( |A8| ) and ( |A12j )) is not literally valid, because each image term in ( |A12b| ) is exactly 



as large in modulus as the "main" term! 

(ii) In what sense does such an expansion correspond to an asymptotic expansion for 
E\(x,y) as A — ► +oo? Formulas (0) would follow immediately from (P if 

oo 

E x (x,y)~\ d / 2 Y,u n \- n/2 (3) 

n=0 

with a n an appropriate multiple of a n . The converse implication from (|[) to (|3|), however, 
is generally not valid beyond the first ( "Weyl" ) term. (For example, in (|Alla| ) or any other 



2 



discrete eigenvector expansion the E\ is a step function; its growth is described by «o but 
there is an immediate contradiction with the form of the higher terms in @.) It has been 
known at least since the work of Brownelli'i that ([3]) is, nevertheless, correct if somehow 
"averaged" over sufficiently large intervals of the variable A. That is, it is valid in a certain 
distributional sense. Hormander&i reformulated this principle in terms of literal asymptotic 
expansions up to some nontrivial finite order for each of the Riesz means of E\ . Riesz means 
generalize to (Stieltjes) integrals the Cesaro sums used to create or accelerate convergence 
for infinite sequences and series (see Section [II]). 

(iii) // an ordinary asymptotic expansion for G does not exist, does an expansion exist in 
some "averaged" sense? We noted above that the Schwinger-DeWitt expansion 

oo 

U(t,x,y) ~ (Amt)- d ^e l ^ 2 / At £ a n (x, y)(it) n ' 2 (4) 

n=0 

is not a true asymptotic expansion under the most general conditions. Nevertheless, this 
expansion gives correct information for the purposes for which it is used by (competent) 
physicists. Clearly, the proper response in such a situation is not to reject the expansion as 
false or nonrigorous, but to define a sense (or more than one) in which it is true. At this point 
we cannot go into the uses made of the Schwinger-DeWitt expansion in renormalization in 
quantum field theory (where, actually, if is a hyperbolic operator instead of elliptic). We 
can note, however, that if U is to satisfy the initial condition in (|A2|) , then as t J, the 
main term in ( |A12b| ), which coincides with the whole of ( |A8b| ), must "approach a delta 



function", while the remaining terms of (|A12b ) must effectively vanish in the context of 



the integral lim t ^ Jo U(t, x, y)f(y) dy. These things happen by virtue of the increasingly 
rapid oscillations of the terms, integrated against the fixed test function f(y). That is, this 
instance of (|j) is literally true when interpreted as a relation among distributions (in the 
variable y). All this is, of course, well known, but our purpose here is to examine it in a 
more general context. We shall show that the situation for expansions like (§) is much like 
that for (||): They can be rigorously established in a Riesz-Cesaro sense, or, equivalently, in 
the sense of distributions in the variable t. This leaves open the next question. 

(iv) If an asymptotic expansion does not exist pointwise, does it exist distributionally in 
x and/or y; and does the spectral expansion converge in this distributional sense when it 
does not converge classically? What is the connection between this distributional behavior 
and that in t? Such formulas as ( |A8a| ), ( |A10| ), ( |A12a|) , (|A14a| ) are not convergent, but 



only summable or, at most, conditionally convergent. The Riesz-Cesaro theory handles the 
summability issue, and, as remarked, can be rephrased in terms of distributional behavior 
in t. However, one suspects that such integrals or sums should be literally convergent in the 
topology of distributions on M. or M. x M.. 

This interpretation is especially appealing in the case of the Wightman function (see 
(|A"4|)-([A~5D, ( |Aiq ), ( |A14| )). To calculate observable quantities such as energy density in 



quantum field theory, one expects to subtract from W(t,x, y) the leading, singular terms 
in the limit y — > x; those terms are "local" or "universal", like the a n in the heat kernel. 
The remainder will be nonlocal but finite; it contains the information about physical effects 
caused by boundary and initial conditions on the field. (See, for instance, Ref. [| Chapters 
5 and 9.) The fact that this renormalized W(t,x,x) is finite does not guarantee that a 
spectral integral or sum for it will be absolutely convergent. Technically, this problem may 
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be handled by Riesz means or some other definition of summability; but in view of the 
formulation of quantum field theory in terms of operator-valued distributions, one expects 
that such summability should be equivalent to distributional convergence on M.. It was, in 
fact, this problem that originally motivated the present work and a companion paper .0. 

A fully satisfactory treatment of these issues cannot be limited to the interior of Ai; 
it should take into account the special phenomena that occur at the boundary. These 
questions are related to the "heat content asymptotics" recently studied by Gilkey et al.BllI 
and McAvityMB (A longer reference list, especially of earlier work by Van den Berg, is 
given by Gilkey in Ref. |T5| .) 

(v) Is the expansion "local" or "global" in its dependence on H? We have already 
encountered this issue in connection with the Wightman function, but it is more easily 
demonstrated by what we call the "cylinder kernel" T(t, x, y), defined by ( |A3| ) . Examination 
of (|A9b| ) and (| A 1 3b|) (| A 1 3 cj) shows that T has a nontrivial power-series expansion in t, 



which is different for the two cases (A4 = R and (0, 7r)). (See Ref. ID for more detailed 



discussion.) More generally speaking, T(t, x, x) differs in an essential way from K(t, x, x) in 
that its asymptotic expansion as t J. is not uniquely determined by the coefficient functions 
(symbol) of H, evaluated at x. T(t, x, x) can depend upon boundary conditions, existence 
of closed classical paths (geodesies or bicharacteristics), and other global structure of the 
problem. In terms of an inverse spectral problem, the asymptotic expansion of T gives more 
information about the spectrum of H and about E\(x,y) than that of K does. (Of course, 
the exact heat kernel contains, in principle, all the information, as it is the Laplace transform 
of E\ .) We shall investigate the issue of locality for a general Green function ([!]). 

In summary, the four basic examples introduced in the Appendix demonstrate all pos- 
sible combinations of pointwise or distributional asymptotic expansions with local or global 
dependence on the symbol of the operator: 





Pointwise 


Distributional 


Local 


Heat 


Schrodinger 


Global 


Cylinder 


Wightman 



(5) 

In this paper we show that the answers to questions (i) and (iii), and the distinction 
between the columns of the table above, are determined by the behavior of g at infinity: If 

g^){t) = 0(f"™) ast^+oo for some 7 eR (6) 

(i.e., g has at infinity the behavior characterizing the test-function space K, — see Sec. ||), 
then the answer to (i) is Yes. On the other hand, when g is of slow growth at infinity 
but does not necessarily belong to /C, then the expansion holds in the distributional sense 
mentioned in (iii). 

The answer to (v), and the distinction between the rows of the table, depend on the 
behavior of g at the origin. If g(t) has an expansion of the form J2^Lo a nt n as t j (even in 
the distributional sense) then the expansion of G(t,x,y) is local. However, if the expansion 
of g(t) contains fractional powers, logarithms, or any other term, then the locality property 
is lost. This subject is treated from a different point of view in Ref. [II]. 

We hope to return to question (iv) in later work. 
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Our basic tool is the study of the distributional behavior of the spectral density e\ = 
dEx/dX of the operator H as A — > oo. We are able to obtain a quite general expansion of e\ 
when H is self-adjoint. Using the results of a previous paper ,@ one knows that distributional 
expansions are equivalent to expansions of Cesaro-Riesz means. Thus our results become 
an extension of those of Hormander.0i They sharpen and complement previous publications 
by one of us.@§§ 

The other major tool we use is an extension of the "moment asymptotic expansion" to 
distributions, as explained in Section |V[ 

The plan of the paper is as follows. In Section [TJ we give some results from Ref. |T^ that 
play a major role in our analysis. In particular we introduce the space of test functions JC 
and its dual /C', the space of distributionally small generalized functions. 

In the third section we consider the distributional asymptotic expansion of spectral de- 
compositions and of spectral densities. Many of our results hold for general self-adjoint 
operators on a Hilbert space, and we give them in that context. We then specialize to the 
case of a pseudodifferential operator acting on a manifold and by exploiting the pseudolo- 
cality of such operators we are able to show that the asymptotic behavior of the spectral 
density of a pseudodifferential operator has a local character in the Cesaro sense. That such 
spectral densities have a local character in "some sense" has been known for years;0H^ here 
we provide a precise meaning to this locality property. 

In the next section we consider two model examples for the asymptotic expansion of 
spectral densities. Because of the local behavior, they are more than examples, since they 
give the asymptotic development of any operator locally equal to one of them. 

In Section |V| we show that the moment asymptotic expansion, which is the basic building 
block in the asymptotic expansion of series and integrals,^ can be generalized to distribu- 
tions, giving expansions that hold in an "averaged" or distributional sense explaining, for 
instance, the small-t behavior of the Schrodinger propagator. 

In the last two sections we apply our machinery to the study of the asymptotic expansion 
of general Green kernels. In Section [VI] we show that the small-t expansion of a propagator 
g(tH) that corresponds to a function g that has a Taylor-type expansion at the origin is 
local and that it is an ordinary or an averaged expansion depending on the behavior of g 
at infinity: If g G JC then the regular moment asymptotic expansion applies, while if g ^ /C 
then the "averaged" results of Section [V] apply. In the last section we consider the case when 
g does not have a Taylor expansion at the origin and show that in that case g(tH) has a 
global expansion, which depends on such information as boundary conditions. 

Some applications of both of the main themes of this paper have been made elsewhere,^ 
most notably a mathematical sharpening of the work of Chamseddine and Corniest on a 
"universal bosonic functional". 



II. PRELIMINARIES 

The principal tool for our study of the behavior of spectral functions and of the associated 
Green kernels is the distributional theory of asymptotic expansions, as developed by several 
authors.&B The main idea is that one may obtain the "average" behavior of a function, in 
the Riesz or Cesaro sense, by studying its parametric or distributional behavior.0 
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In this section we give a summary of these results. We also set the notation for the 
spaces of distributions and test functions used. 

If A4 is a smooth manifold, then T>(Ai) is the space of compactly supported smooth 
functions on M., equipped with the standard Schwartz topology.cJca Its dual, T>'(A4), is the 
space of standard distributions on A4. The space £(Ai) is the space of all smooth functions 
on M., endowed with the topology of uniform convergence of all derivatives on compacts. 
Its dual, £'{M), can be identified with the subspace of V\M) formed by the compactly 
supported distributions. Naturally the two constructions coincide if M. is compact. 

The space iS^R") consists of the tempered distributions on R n . It is the dual of the 
space of rapidly decreasing smooth functions iS(R n ); a smooth function (f> belongs to 5(R ?1 ) 
if D a <fi(x) = o(\x\~°°) as \x\ — > oo, for each a G N n . Here we use the usual notation, 
jja _ Q\ a \ jQx® 1 ■ ■ ■ dx% n , \oc\ — oi\ + • — h a n ; o(x~°°) means a quantity that is o(x _/3 ) for all 
/?GR. 

A not so well known pair of spaces that plays a fundamental role in our analysis is /C(R n ) 



and JC'(R n ). The space /C was introduced in Ref. |29|. A smooth function belongs to K, q 
if D a (j)(x) = 0(\x\ q ~\ a \) as \x\ — > oo for each a G N n . The space K. is the inductive limit of 
the spaces JC q as q — » oo. 

Any distribution / G /C'(R) satisfies the moment asymptotic expansion, 

/(Ax) ~ ■!?,■+! asA^oo, (7) 

3=0 J- A 

where fij = {f{x),x^) are the moments of /. The interpretation of (0) is in the topology 
of the space /C; observe, however, that there is an equivalence between weak and strong 
convergence of one-parameter limits in spaces of distributions, such as K! . 

The moment asymptotic expansion does not hold for general distributions of the spaces 
V or S' . Actually, it was shown recentlyS that any distribution / G T>' that satisfies the 
moment expansion (|7|) for some sequence of constants {/i, } must belong to K' (and then the 
j are the moments). 

There is still another characterization of the elements of K, 1 . They are precisely the 
distributions of rapid decay at infinity in the Cesaro sense. That is why the elements of /C' 
are referred to as distributionally small. 

The notions of Cesaro summability of series and integrals are well known.S In Ref. [16 



this theory is generalized to general distributions. The generalization includes the classical 
notions as particular cases, since the behavior of a sequence {a n } as n —>■ oo can be studied 
by studying the generalized function J2^=o a n<5(x — n )- The basic concept is that of the 
order symbols in the Cesaro sense: Let / G £>'(R) and let (3 G R\ {—1, —2, —3, . . .}; we say 
that 

f(x) = 0{x p ) (C) as x -> oo, (8) 

if there exists iV G N, a function F whose Ath derivative is /, and a polynomial p of degree 
N — 1 such that F is locally integrable for x large and the ordinary relation 

F{x) = p(x) + 0(x p+N ) as x -> oo (9) 

holds. The relation f(x) = o(ar) (C) is defined similarly by replacing the big O by the 
little o in (M. 
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Limits and evaluations can be handled by using the order relations. In particular, 
lim^oo f(x) = L (C) means that f(x) = L + o(l) (C) as x — > oo. If / e V' has 
support bounded on the left and G £, then in general the evaluation (f(x),(f)(x)) does 
not exist, but we say that it has the value S in the Cesaro sense if lim^oo G(x) = S (C), 
where G is the primitive of f<ft with support bounded on the left. The Cesaro interpretation 
of evaluations (f(x), 4>(x)) with supp/ bounded on the right is similar, while the general 
case can be considered by writing / = fx + f 2 , with supp fi bounded on the left and supp f% 
bounded on the right. 

The main result that allows one to obtain the Cesaro behavior from the parametric 
behavior is the following. 

Theorem 2.1. Let f be in V with support bounded on the left. If a > —1, then 

f(x) = 0(x a ) (C) as x -> oo (10) 

if and only if 

/(Ax) = 0{X a ) as A -> oo (11) 

distributionally. 

When — {k + 1) > a > —(k + 2) for some k G IN, holds if and only if there are 
constants /i , . . . , Hk su °h that 

/N^t^ + W (12) 



j=0 



j\ A-?' +1 



distributionally as A — > oo. 
Proof: See Ref. 0. ■ 

The fact that the distributions that satisfy the moment asymptotic expansion are exactly 
those that satisfy f(x) = 0(x~°°) (C) follows from the theorem by letting a —>■ — oo. Thus 
the elements of JC' are the distributions of rapid distributional decay at infinity in the Cesaro 
sense. Hence the space K! is a distributional analogue of S. We apply this idea in Section [V], 
where we build a duality between S' and K! . 

Another important corollary of the theorem is the fact that one can relate the (C) 
expansion of a generalized function and its parametric expansion in a simple fashion. Namely, 
if {dj} is a sequence with 3?e ctj \ — oo, then 



f(x) ~ a 3 xa ' (°) as x ^ oo (13) 

j=0 



if and only if 



00 00 (-i)j u .xU)( T ) 

/(Ax) ~ J2 + E ^ ( 14 ) 

i=o i=o J- A 

as A — >• oo, where the fij are the (generalized) moments of / and where 

g a (x)=xl if -1,-2,-3,..., (15) 
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while in the exceptional cases g a is a finite-part distribution:^ 

g~k(x) = V.f. (x(x)x~ k ) if k = 1,2,3,..., (16) 

X being the Heaviside function, the characteristic function of the interval (0, oo). Notice 
that 

g a (Xx) = X a g a (x), a ^ -1,-2,-3,..., (17) 

9-k{)*) = ^r L + ± \ k _ 1)lxk 1 J > fe = l,2,3,.... (18) 



III. THE ASYMPTOTIC EXPANSION OF SPECTRAL DECOMPOSITIONS 

Let 7i be a Hilbert space and let if be a self-adjoint operator on 7f , with domain X . 
Then if admits a spectral decomposition {^aIaL-oo • The {-Ea} is an increasing family of 
projectors that satisfy 

/ = / dE\ , (19) 
where i is the identity operator, and 

H = I XdE x (20) 

in the weak sense, that is, 

(Hx\y) = I Xd(E x x\y), (21) 



oo 



oo 



OO 



for x £ X and y G where is the inner product in 7i. 

Perhaps more natural than the spectral function E\ is the spectral density e x = dE\/dX. 
This spectral density does not have a pointwise value for all A. Rather, it should be un- 
derstood as an operator-valued distribution, an element of the space V(R, L(X, TO)). Thus 
(H)-© become 

i=(e A ,l) (22) 



H=(e x ,\), (23) 

where (/(A), 0(A)) is the evaluation of a distribution /(A) on a test function 0(A). 

The spectral density e x can be used to build a functional calculus for the operator H. 
Indeed, if g is continuous and with compact support in R then we can define the operator 
g(H) e L(X,H) (extendible to L(H,H)) by 

g(H) = (e x ,g(X)). (24) 
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One does not need to assume g of compact support in (|24D , but in a contrary case the domain 
of g(H) is not X but the subspace M g consisting of the x G Ti for which the improper integral 
((e\x\y), g(\)) converges for all y ETC. 

One can even define f(H) when / is a distribution such that the evaluation (e\, /(A)) is 
defined. For instance, if E\ is continuous at A = Ao then E\ = x(Ao — H) where x is again 
the Heaviside function. Differentiation yields the useful symbolic formula 

e x = 6(X-H). (25) 

oo 

Let X n be the domain of H n and let X^ = P| X n . Then 

n=l 

(e x ,X n ) = H n (26) 

in the space L(X 00 ,H). But, as shown recently!!! a distribution / G T>'(R) whose moments 
(f(x),x n ), n G N, all exist belongs to /C'(R), that is, is distributionally small. Hence, e\ , as 
a function of A, belongs to the space /C'(R, L(A' 00 , H)). Therefore, the asymptotic behavior 
of e\ a , as cr — > oo, can be obtained by using the moment asymptotic expansion: 



oo 



(-l)"^(")(Aa) 

e A<x ~ 2^ j as (x ^ oo, (27) 



„=o n! 



while eA vanishes to infinite order at infinity in the Cesaro sense: 

e A = o(|A|-°°) (C) as |A| -> oo. (28) 

The asymptotic behavior of the spectral function is obtained by integration of ( |27|) 
and by recalling that lim £\ = 0, lim E\ = I. We obtain 

A^— oo A— *oo 

oo / i \n+l TTn+1 s(n) ( \ _\ 

E\ ~ x(Acr)/ + ^ lj d W asa^oo. (29) 

n=0 l n + X J- 

Similarly, the Cesaro behavior is given by 

E x = I + o(\-°°) (C) asA^oo, (30) 



E\ = o(|A|-°°) (C) as A -> -oo. (31) 

These formulas are most useful when if is an unbounded operator. Indeed, if H is 
bounded, with domain X — 7i, then e\ = for A > and -E^ = for A < — \\H\\, E\ = I 
for A > \\H\\, so (|2"B|), (PPf), and ( PT|) are trivial in that case. 

In the present study we are mostly interested in the case when H is an elliptic differential 
operator with smooth coefficients defined on a smooth manifold M.. Usually TC = L 2 (M) 
and X is the domain corresponding to the introduction of suitable boundary conditions. 
Usually the operator H will be positive, but at present we shall just assume H to be self- 
adjoint. 
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In this case the space V(M) of test functions on M. is a subspace of X^ . Observe also 
that the operators K acting on T>(M) can be realized as distributional kernels k(x,y) of 
V'(M x M) by 

{K<P){x) = (k{x iy )^{y)) y . (32) 

In particular, 5(x — y) is the kernel corresponding to the identity /, and H5(x — y) is the 
kernel of H. The spectral density e\ also has an associated kernel e(x,y; A), an element of 
V(1R, V(Ai x M)). Since H is elliptic it follows that e(x, y; A) is smooth in (x, y). Warning: 
Much of the literature uses u e(x,y; A)" for what we call ll E(x,y; A)". 

The expansions (|T7|)-(|31~|) will hold in X^ and thus, a fortiori, in T>(M). Hence 

e(x,y; Xa) ~ 2^ j as a ^ oo, (33) 

n=o n - 

£;(x,i/;A(t)~x(A(t)5(x-i/) + ^^ ) y; ^ asa^oo, (34) 

n=o (71 + 1J. 

in the space V x (R,V xy (M x A4)). Furthermore, 

e(x,y;A) = o(|A| _0 °) (C) as |A| -> oo, (35) 



E(x,y;X) = 5(x - y) + o^- 00 ) (C) as A -> oo, (36a) 

E(x,y; A) = o(|A|~°°) (C) as A -> -oo, (36b) 

in the space £>'(.M xM). 

Actually, an easy argument shows that the expansions also hold distributionally in one 
variable and pointwise in the other. (For instance, fl35|) says that if y is fixed and G T>(M) 
then (e(x,y; \),<p(x)) = o(|A|-°°) (C) as |A| -> oo.) 

That the expansions cannot hold pointwise in both variables x and y should be clear since 
we cannot set x — y in the distribution 5(x — y). And indeed, e(x, x; A) is not distributionally 
small. However, as we now show, the expansions are valid pointwise outside of the diagonal 
of M x M. 

Indeed, let U, V be open sets with U n V = 0. If / G D'(.M) and G £>(R), 
then is a smoothing pseudodifferential operator, so <$>(H)f is smooth in .M. Thus, 

(e{x,y;\),f{x)g(y)<j>{\)) = (<j>(H)f{x),g(x)) is well-defined if / G D'(.M), supp/ CU,ge 
V'(M),suppg C 1/. Therefore e(a:,y; A) belongs to Z>'(R,£(C7 x 1/)). But 

(e(x, y; A), f(x)g(y)X n ) = (H n f(x),g(x)) = 0, (37) 

thus e(x, y;X) actually belongs to /C'(R, £{U x V)); that is, it is a distributionally small 
distribution in that space whose moments vanish. Therefore 

e(x, y; Xa) = o(a~°°) as a — > oo, (38) 
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E(x,y; Xa) = x(Xa)5(x — y) + o(a °°) as a — > oo, (39) 



in the space /C'(R, £{U x V)). Similarly, (|35|) — (|36|) also hold in £(U x V). Convergence in 
£([/xl^) implies pointwise convergence on UxV, but it gives more; namely, it gives uniform 
convergence of all derivatives on compacts. Thus (pop, (|36|), (|38|), and ( |39"D hold uniformly 
on compacts of U x V and the expansion can be differentiated as many times as we please 
with respect to x or y. 

Example. Let Hy = —y" considered on the domain X = {y G C 2 [0, tt] : y(0) = y{ji) = 
0} in L 2 [0, 7r]. The eigenvalues are X n = n 2 , n = 1, 2, 3, . . ., with normalized eigenf unctions 
4>n(x) — \ sin no; . Therefore, 

2 00 

e(x, ?/; A) = — y2 sin nx sin ny 5(X — n 2 ) , (40) 

n n=l 

where < x < 7r, < y < ti. Then 

2 ~ . ™6<M(x-y)6M(\) 

— y j sin nx sm nyo(Xa — n ) ~ y j q — — as cr — > oo (41) 

in V'(n,V'{{0,ir) x (0,tt))), while 
2 00 

— V] sinnx sinm/<5(Acr — n 2 ) = o(cr~°°) as o — > oo (42) 

^ n=l 

if x and ?/ are fixed, x ^ y. On the other hand, 

2 00 

e(x, x; A) = — sin 2 nx S(X — n 2 ), (43) 

thus if < x < 7r, 



e(x, x; Aa) = — — cos2nx) S(Xa — n 2 ^ 



n=l 



1 00 1 

- V 5(Acr - n 2 ) H 5(X) + o(cr~°°) as cr -> oo, 

7T ^ 27TCT 



because the generalized function Y^=\ cos 2nx S(X— n 2 Hs distributionally small if < x < 7r, 
with moments /xq = —1/2 and /j& = for > 1, sincetll 



OO 

^]cos2nx = -- (C), 

n=l 1 



^n 2fc cos2nx = (C), fc = 1, 2, 3, . . . . 

n=l 

But (Ref. |6|, Chapter 5) 

00 1 r°° 1 
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as e — > + if G S, thus 

X —1/2 

e(x, x; Act) = - — A + + o(cr _ °°) as cr — > oo. (45) 
It is then clear that e(x, x; A) is not distributionally small; rather, 

e(x,x;X) = - /2 + o(\-°°) (C) as A - oo, (46) 

Z7rA 

that is, e(x,x; A) ~ (l/27r)A -1 / 2 , as A — > oo, in the Cesaro sense. ■ 

Neither is the spectral density e(x, y; A) distributionally small at the boundaries, as 
follows from the heat content asymptotics of Refs. 0,^. That there is a sharp change of 
behavior at the boundary can be seen from the behavior of the spectral density e(x, x; A) 
given by (p|). Indeed, if < x < n then e(x, x; A) = (l/2n)\~ 1/2 + o(A _ °°) (C), but when 
x = or x = Ti then e(0, 0; A) = e(n, n; A) = 0. 

It is important to observe that in the Cesaro or distributional sense, the behavior at 
infinity of the spectral density e(x, y\ A) depends only on the local behavior of the coefficients 
of H. That is, if H\ and H2 are two operators that coincide on the open subset U of Ai and 
if ei(x, y; A) and e2(x,y; A) are the corresponding spectral densities, then 

ei(x, y; crA) = e%(x, y; aX) + o(a -00 ) as a — > 00 (47) 

in D'(U x U). This follows immediately from (|33|). In fact, it follows from Theorem 7.2 that 

ei(x,y; A) = e 2 (x, y; A) + o(A"°°) (C) as A 00, (48) 

pointwise on (x,y) £ U xU (even on the diagonal!). More than that, (fS|) holds in the space 
£{U x U), so that it is uniform on compacts of {7. These results are useful in connection 
with the suggestionElhii to replace a general second-order operator H by another, H , that 
agrees locally with H and for which the spectral density can be calculated. In the next 
section we treat two special classes of operators where this idea has been implemented. 

Example. The spectral density for the operator — y" on the whole real line is 

/ u x(A)cosA 1/2 (x-?/) 
e 1 {x,y;X) = ' ( 49 ) 



as can be seen from ( A7a ) and ( |A8a| ). Therefore, comparison with (HO) yields 



2 00 cos \ l l 2 (r — v) 

-Ysmnxsmny8(\-n 2 ) = U> + ofA" 00 ) (C) as A -> 00. (50) 

7T ^ 2ttA 1 /2 



In particular, if we set x — y we recover 

Formula fl50"D is uniform in compacts of (0, n) x (0, 7r) but ceases to hold as x or y 
approaches or 71. For instance, if y = 0, the left side vanishes while [cf. (|33|)1 

X^cos^A) 1 / 2 ^ 8(x)S(X) 8"(x)8'(X) 
2tt(ctA) 1 /2 a a 2 

as a — > 00. ■ 
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IV. SPECIAL CASES 



In this section we give two model cases for the asymptotic expansion of spectral densities. 
They are not just examples, since according to the results of the previous section, the spectral 
density of any operator locally equal to such a model case will have the same behavior at 
infinity in the Cesaro sense. 

Let us start with a constant-coefficient elliptic operator H defined on the whole space 
R n . Then H admits a unique self-adjoint extension (which we also denote as H), given as 
follows. Let p = o~(H) be the symbol of H (i.e., H = p(—id)). Then the spectral function is 
given by 



E(x,y;X) = —- e*W d£, (51) 

(2tt)" y P (o<A 



so that the spectral density can be written as 



e(x, vA) = j^y (e l{x ~ yH , *(P(0 - A)) • (52) 



For the definition of 5(f(x)) see Refs. 

To obtain the behavior of e(x, y; A) as A — > oo in the Cesaro or in the distributional 
sense, we should consider the parametric behavior of e(x, y; aX) as o — > oo. Setting e = 1/er 
and evaluating at a test function 0(A), one is led to the function 

$(e) = < e (x,i/;A) s 0(eA)) A . (53) 

But in view of Q52D we obtain 

m = j^(e l(x ' yH ,<P(ep(0)) r (54) 

When x ^ y are fixed, e l ^ x ~ y ''^ is distributionally small as a function of £. This also holds 
distributionally in (x,y). Thus the expansion of ([54]) follows from the following lemma. 

Lemma 4.1. Let f G K.'(M n ), so that it satisfies the moment asymptotic expansion 

{-l)\ k \^ k DH{x) 

/(Ax)~ Y. ju x \ k \ +n «sA^oo, (55) 

where fik = (f{x),x k ) , k e IN* 1 , are the moments. Then if p is an elliptic polynomial and 
0G/C, 

iff \M f K\ ^ (f(x),p(x) m ) 0M(O)g m 

(/(x), 0(ep(a;))) ~ 2^ j as e 0. (56) 

n=0 m - 

Proof: The proof consists in showing that the Taylor expansion in e, 

^)) ^ m){0) f rem + O(e^), (57) 
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not only holds pointwise but actually holds in the topology of /C(R n ). But the remainder 
in this Taylor approximation is 

^(6p(x))p(x) N +^ 
K N {x, e) (JV + 1)! ' ^ ^ 

for some 6 £ (0, e). Since G /C, there exists q e R such that 4>^\x) = 0{\x\ q ^^) as x — > oo. 
If p has degree m it follows that 

|i?AT(x,£)i< (Vil)! (59) 

for some constant M, and the convergence of the Taylor expansion in the topology of the 
space fC follows. ■ 

Thus, applying (|||) with f(x) = e l ( x ~ y ^ for x ^ y or distributionally in (x, y), we obtain 

_J_ ~ (e^)^,p(Q fc )0( fc )(O) g fc 

" (27T)" fr k\ 



or 



^^■^fc^W (60) 

fc=o k - 



Therefore, 



« (-l) k H k 6(x-y)S^(X) 
e{x,y; Act) ~ ^ as a ^ oo, (61) 

in accordance with the general result. 

Observe also that if Hi is any operator corresponding to the same differential expression, 
considered in some open set Ai with some boundary conditions, then its spectral density 
ei(x, y; A) satisfies 

ei (x,y;A) = -^(e^H <y(p(£) - A)) + o(A°°) (C) as A - oo. (62) 

Example. Let M. be a region in R ra and let if be any self-adjoint extension of the 
negative Laplacian —A obtained by imposing suitable boundary conditions on M.. Let 
e M{x, Wi X) be the spectral density. Then 

e M (x, y; A) = J^yMltf ~ A), e^*) + (A"~) (C). (63) 

We now use the one-variable formula 

8(x — x ) 



*(/(*)) 



\fU)\ 
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valid if / has a single zero at xq , and pass to polar coordinates £ = ru, where r = |£|, 
uj = (ui, . . . ,u n ) satisfies |u;| = 1, and = r n ~ l drda(u>), to obtain 

(5(|£| 2 - A), e^-wK) = -f— / / (5{r 2 - A), e ^-v)^-i dr da (^) 



(27r) n 1 n (27r) n J|o.]=i Jo 



i\ 1 ' 2 (x— y)-u 



2(2vr)™ 7h=i 



da{ui) 



iX^u> 1 \x-y\ d / n 

2(27r) n 7| w |=i V ; 

\n/2-l o_(n-l)/2 -i 

2(27r) n r(2fi) 7-1 V ; 

^- 1/2 Jn,2-l{\ 1/2 \x-y\) 
2n/2+l 7J -n/2| a; _ ^|n/2-l ' 

where J p (x) is the Bessel function of order p. Therefore 

e M (x,y;X) = 2n/2+1 ^ x _ y ^ 2 _ x + o(X ) (C) as A - oo, (64) 

uniformly over compacts of M. x M.. ■ 

Our second model is an ordinary differential operator H with variable coefficients, as 
treated in Refs. 18-20. There are two major simplifications in this one-dimensional case. 
First, the Weyl-Titchmarsh-Kodaira theoryBi! expresses the spectral density as 

i 

e(x,y;\)d\ = ^ 4>\j(x) dp? k {\) if>xk(v) , (65) 

j,k=0 

where d^ k are certain Stieltjes measures supported on the spectrum of H, and ip\j are the 
classical solutions of Hip — Xip with the basic data 

^xo(x ) = 1, if>\o(x ) = 0, 

4>\i(x ) = 0, ip' xl {xo) = 1, (66) 

at some x G M.. Thus 

BE 

V 00 (\) = E{x ,x ;\), /i 01 (A) = _(a;o,xo;A), 

dy 

dE d 2 E 

yU 10 (A) = -fofro, x ; A), /i n (A) = q^q-( x o, x ; A). (67) 

Second, the eigenfunctions can be approximated for large A quite explicitly by the 
phase-integral (WKB) method. (Thirdly, but less essentially, there is no loss of generality 
in considering 

H=-^ + V(x), (68) 
15 



since the general second-order operator can be reduced to this form by change of variables.) 

In Ref. [18] the phase-integral representation of the eigenfunctions was used to obtain in 
a direct and elementary way the expansion 

1 oo 

d/i jk (\) ~ - V p> k (x ) u 25 ^- 2n duj, (69) 

where A = uj 2 and 

p8° = i, pT = \v, pI° = 1(-v" + 3v 2 ), 

^ = 1, Pi=~\ V > pl^liV'-SV 2 ), (70) 

p 10 = p 01 = -— (p 00 ). 

Formula (|69| ) is a rigorous asymptotic expansion when M. = R and is a C°° function 
of compact support. The relevance of ( |69"D in more general cases, where it is certainly not 



a literal pointwise asymptotic expansion, was discussed at length in Ref. [18]; the results of 
the present paper simplify and sharpen that discussion by showing that the error in ( jff^J 
is 0(A~°°) in the (C) sense for any operator locally equivalent to one for which (j^U ) holds 
pointwise. 

V. POINTWISE AND AVERAGE EXPANSIONS 

Let / G /C(R). Since the elements of /C'(R) are precisely the distributionally small 
generalized functions, it follows that / satisfies the moment asymptotic expansion; that is, 

/(Ax) ~ ^ as A -> oo, (71) 



where 



/' ; (/(.'')•.'■')• jeN, (72) 



are the moments. 

The moment asymptotic expansion allows us to obtain the small-t behavior of functions 
G(t) that can be written as 

G(t) = (f(x),g(tx)), (73) 
as long as g G K. Indeed, (|7T|) gives 

G(«) = f a si ^0. (74) 
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Naturally, this would be valuable if /(A) = e(x, y; A) is the spectral density of the elliptic 
differential operator H and G(t,x,y) = (e(x,y; A), g(Xt)) is an associated Green kernel. 

However, we emphasize that the derivation of fl74|) holds only when g G fC. What if 
g /C ? A particularly interesting example is the kernel U(t,x,y) = (e(x,y;X), e~ lXt ) that 
solves the Schrodinger equation 

dU 

i— = HU, t>0 (75a) 
dt ( J 

with initial condition 

U(0+,x,y) = 6{x-y). (75b) 
In this case g(x) = e~ tx is smooth, but because of its behavior at infinity, it does not belong 



to JC. We pointed out in the introduction, however, that (7j) is still valid in some "averaged" 
sense. 

Indeed, we shall now show that formula (f73|) permits one to define G(t) as a distribution 
when instead of asking g G JC we assume g to be a tempered distribution of the space S' 
which has a distributional expansion at the origin. We then show that (|7i|) holds in an 
averaged or distributional sense. The fact that the space of smooth functions K, is replaced 
by the space of tempered distributions is not casual: the distributions of S' are exactly 
those that have the behavior at oo of the elements of K, in the Cesaro or distributional 
sense. Indeed, we have 

Lemma 5.1 Let g G S'(M). Then there exists a G M such that 

g (n \\x) =0(\ a - n ) as A^oo, (76) 

distributionally. 

Proof: See Ref. |l^, where it is shown that (|76|) is actually a characterization of the 
tempered distributions. ■ 

Let g G (S'(R) and let a be as in (|76|). If G 5(R) then the function $ defined by 

<S>{x) = (g(tx),<t>(t)) (77) 
is smooth in the open set (— oo, 0) U (0, oo) and, because of (fTS]), satisfies 

®( n \x) = 0(|a;r- n ) as |ac| -> oo. (78) 
It follows that we can define G(t) = (f(x),g(tx)) as an element of 5'(R) by 

(G(t),0(t)> = <^)>, (79) 

whenever f E JC' and ^ supp /. 

When G supp / then ( [79]) cannot be used unless $ is smooth at the origin. And in 
order to have $ smooth we need to ask the existence of the distributional values g^(0), 
n = 0,1,2,.... 

Recall that following Lojasiewicz,@ one says that a distribution h G T>' has the value 7 
at the point x = Xq , written as 
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h(x ) =7 in £>', (80) 

if 

lim h(x + ex) = 7 (81) 
distributionally; that is, if for each (ft E P 

/oo 
0(x) dx. (82) 
-00 

It can be shown that Ji(xq) = 7 in T>' if and only if there exists a primitive h n of some 
order n, hffi = h, which is continuous in a neighborhood of x = Xq and satisfies 

7(37 — Xo) n 

h n (x) = ; h o(|a; — a; | n ), as x — > x . (83) 

n! 

In our present case, we need to ask the existence of the distributional values g^ n '(0) = a n 
for fieN. We can then say that g(x) has the small- X "averaged" or distributional expansion 

d( x ) ~ X! n ' as x — > 0, in V', (84) 

n=0 n - 

in the sense that the parametric expansion 

00 a e n x n 

g{ex) ~ 51 — ' as e -> 0, (85) 

n=0 n - 

holds, or, equivalently, that 

(g(ex), cj>{x)) ~ X)-T ( / ^ ( 86 ) 

n=0 n! VJ -°° 7 

for each (ft E T). 

Lemma 5.2. Lei g E S' be such that the distributional values g^ n \0) = a n ,in V , exist 
for n E N. Let <p E S and put $(x) = (g(tx), 4>{t)). Then $ G /C. 

Proof: Indeed, $ is smooth for x 7^ 0, but since the distributional values p^(0) = a n 
exist, it follows that $(x) ~ S^Lo ^n 2 -™ as 1 -> 0, where b n = (a n /n\) f^° 00 x n 4>(x) dx. Thus 
$ is also smooth at x — 0. Finally, let a be as in fl?D; then & n \x) = 0{\x\ a ~ n ) as \x\ -> 00. 
Hence $ G /C. ■ 

Using this lemma we can give the following 

Definition. Let f G /C. Lei g E S' have distributional values g( n '(0), n E IN. Then we 
can define the tempered distribution 

G(t) = (f(x),g(tx)} (87) 

by 

(G{t),<t>{t)) = (f(x),9(x)), (88) 
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where 



$(x) = (g(tx),<j>(t)), (89) 

In general the distribution G(t) is not smooth near the origin, but its distributional 
behavior can be obtained from the moment asymptotic expansion. 

Theorem 5.1. Let f G KJ with moments fi n = (f(x),x n ). Let g G S' have distribu- 
tional values g^(0) for n G IN. Then the tempered distribution G(t) = (f(x),g(tx)) has 
distributional values G^ n \0), n G FN, which are given by G^ n \0) = fj, n g( n \0) ; and G has the 
distributional expansion 

G{t) ~ f; ^ (W) (°)*" , i n v', ast^O. (90) 

n=0 n - 



Proof: Let G S and let $(x) = (g(tx),(f>(t)). Then 

(G(et),0(t)) = (/(x) l $(ex)>, (91) 
and since $( n )(0) = g^ n \0) Jf^ t n cf){t) dt, the moment asymptotic expansion yields 

(G(e*), <f>{t)) ~ £ ^g^igj ( f°° t «0( t ) ^ £ n as £ ^ 0> (92) 

n=0 n - \J-oo J 

and © follows. ■ 

Before we continue, it is worthwhile to give some examples. 

Example. Let g G S' be such that the distributional values g^ n \0) exist for n G IN. 
Since the Fourier transform g(X) can be written as g(X) = X~ 1 (e ix , g(A _1 x)), and since all 
the moments /i n = (e tx ,x n ) vanish, it follows that gie^ 1 ) = 0(e°°) distributionally as e — > 
and thus #(A) = 0(|A|-°°) (G) as |A| -> oo. Therefore <? G /C'. 

Conversely, if / G /C', then its Fourier transform f(t) is equal to = (f(x),e' ltx ) 
for t ^ 0. Thus /(t) = F(i) + X)? = o % <^') (t) f° r some constants ao , . . . , a n . But the 
distributional values F^(0) exist for n G IN and are given by F(")(0) = i n (f(x),x n ), and 
hence F E K', and it follows that ao = • • • = a n = 0. In summary, / (n) (0) exists in V for 
each n G IN. 

Therefore, a distribution g G S' is smooth at the origin in the distributional sense (that 
is, the distributional values g^ n \0) exist for n G IN) if and only if its Fourier transform g is 
distributionally small (i.e., g G K'). ■ 

Example. Let £ gC with |£| = 1, £ / 1. Then the distribution f{x) = En=-oo C^(x-n) 
belongs to KJ . All the moments vanish: fi k = Y%?=-oo C nk = (G) for fc = 0, 1, 2, . . .. It 
follows that if g G S' is distributionally smooth at the origin, then 

oo 

J2 Cg(nx) = o(x°°) in V as x -> 0. (93) 

n=— oo 

When £ = 1, D^L_oo <K X — n ) does not belong to KJ but X^_oo <K X — — 1 does. Thus, 
if g G S' is distributionally smooth at the origin and ff^ du is defined, then 
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00 / POO \ 

g(nx)=[ g(u) du) x" 1 + o(x°°) in V as x -> 0. (94) 

Actually, many number-theoretical expansions considered in Ref. |37] and Chapter 5 of Ref. ^ 
will hold in the averaged or distributional sense when applied to distributions. ■ 

Many times, supp / C [0, oo) and one is interested in G(t) = (f(x),g(tx)) for t > 
only. In those cases the values of g(x) for x < are irrelevant and one may assume that 
suppg C [0, oo). Since we need to consider $(x) = (g(tx), 4>{t)) for x > only, we do not 
require the existence of the distributional values «7^(0); instead, we assume the existence 
of the one-sided distributional values g^ n \0 + ) = a n for n G N. This is equivalent to asking 
g(ex) to have the asymptotic development 



00 a n e n x 



g{ex)~J2 " , + as^0+; (95) 



n=0 



n 



that is 



(g(ex), (f>(x)) ~ E ^[ (y o dx ) £ " as £ ^ 0+ (96) 



n=0 

for G <S. We shall use the notation 

( x )^g^! in p' asx ^0+ (97) 

n=0 n - 

in such a case. 

Lemma 5.3. Let g £ S' with suppg C [0, oo) and let g^ n \0 + ) = a n exist in V for 
n G IN. Let <fi G <S and pirf $(x) = (g(tx),<p(t)) for x > 0. T/ien $ admits extensions $ to 
iR TOf/i $ G /C(iR). 

Proof: It suffices to show that $ is smooth up to the origin from the right and that it 
satisfies estimates of the form &^(x) = 0(|a;| a_n ) as \x\ — > oo. But the first statement 
follows because g^ n '(0 + ) exists for all n G IN, while the latter is true because of ([76]) . ■ 

From this lemma it follows that when / G K! , supp / C [0, oo), supp g C [0, oo), and the 
distributional values g^ n \0 + ) exist for n G N, then G(t) = (f(x),g(tx)) can be defined as a 
tempered distribution with support contained in [0, oo) by 

(G(t),0(t)) = (/(x),i>(x)), (98) 

where $ is any extension of = (g(tx), 4>(t)), x > 0, such that $ G AT. 

Theorem 5.2. Lei / G /C' witt supp / C [0, oo) and moments \i n = (f(x), x n ) . Let g G 
S' with supp a C [0, oo) have distributional one-sided values g( n \0 + ) for n G IV. Then the 
tempered distribution G(t) = (f(t),g(tx)) defined by ft9l\ ) has distributional one-sided values 
G ( - n \0 + ), n G IN, which are given by G^ n ^(0 + ) = /i„g( n )(0 + ), and G has the distributional 
expansion 

00 ii n^(() + )t n 

g(*)~E I mV ' ast^0 + . (99) 

n=0 n ' 



Proof: Quite similar to the proof of Theorem 5.1 . 
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VI. EXPANSION OF GREEN KERNELS I: LOCAL EXPANSIONS 



In this section we shall consider the small-t behavior of Green kernels of the type 
G(t;x,y) = (e(x,y; X), g(Xt)) for some g G S'. Here e(x,y;X) is the spectral density ker- 
nel corresponding to a positive elliptic operator H that acts on the smooth manifold M.. 

Our results can be formulated in a general framework. So, let H be a. positive self-adjoint 
operator on the domain X of the Hilbert space 7i. Let X^ be the common domain of H n , 
n G IN, and let t\ be the associated spectral density Let g G S' with supp g C [0, oo) such 
that the one-sided distributional values a n = 5 f< - n - ) (0 + ) exist for n G IN. Then we can define 

G(t) = g{tH), t > 0, (100) 

that is, 

G(t) = (e x ,g(tX)), t>0. (101) 

Thus G can be considered an operator- valued distribution in the space iS'(R, L(X OQ , TC)). The 
behavior of G(t) as t — > + can be obtained from the moment asymptotic expansion (j27D 
for e\ . The expansion of G(t) as t —>■ + will be a distributional or "averaged" expansion, 
in general, but when g has the behavior of the elements of JC at oo it becomes a pointwise 
expansion. In particular, if g is smooth in [0, oo), the expansion is pointwise or not depending 
on the behavior of g at infinity. 

Theorem 6.1. Let H be a positive self- adjoint operator on the domain X of the Hilbert 
space TC. Let X^ be the intersection of the domains of H n for n G IN. Let g G S' with 
supp g C [0, oo) be such that the distributional one-sided values 
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W(0+) = a n inV (102) 



exist for n G IN. Let G(t) = g(tH), an element of S' '(M, L(X 00 ,Tt)) with support contained 
in [0, oo). Then G(t) admits the distributional expansion in L{X OQ: l-C), 

oo H n t n 

G (t) ~ E — > as * 0+ > m V '> ( 103 ) 

n=0 n - 

so that the distributional one-sided values G^(0 + ) exist and are given by 

G {n \0 + ) = a n H n mV. (104) 



When g admits an extension that belongs to 1C, (\103j) is an ordinary pointwise expansion 
while the G ( - n \0 + ) exist as ordinary one-sided values. 

Proof: Follows immediately from Theorem 5.2. ■ 

When if is a positive elliptic differential operator acting on the manifold Ai, then The- 
orem 6.1 gives the small-t expansion of Green kernels. Let e(x,y; A) be the spectral density 
kernel and let 

G{t,x,y) = (e{x,y;X),g(tX)), t > 0, (105) 
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be the Green function kernel corresponding to the operator G(t) = g(tH). Then G belongs 
to S'(R)®V(M. x A4), has spectrum in [0, oo), and as t — > + admits the distributional 
expansion 

x, y) ~ X] — — n — ' ast ^° » inp > ( 106 ) 

n=0 n ' 

that is, 

G(et,x,y) ~ 2^ j — as £ — > , (107) 

n=o n - 

d n 

in D'(.M x M). Also, the distributional one-sided values -^G(0 + , x, y) exist for n G IN and 
are given by 

—G(0 + ,x,y) = a n H n 5{x-y) in V. (108) 

If g admits extension to /C, then ( |106| ) and ( |108| ) are valid in the ordinary pointwise sense 
with respect to t (and distributionally in (x,y)). 

Pointwise expansions in (x, y) follow when x y. Indeed, if U and V are open subsets 
of M with U H V = 0, then G belongs to S'(R)<g)£{U x V) and as t -> 0+ we have the 
distributional expansion 

G(t, x, y) = o(t°°), in Z>', as i -> + , (109) 

in x V), and in particular pointwise on x G Z7 and y G V. The expansion becomes 
pointwise in i when g admits an extension to /C. 

These expansions depend only on the local behavior of the differential operator. Let 
Hi and H 2 be two differential operators that coincide on the open subset U of M.. Let 
e±(x, y,X), e2(x,y,X) be the corresponding spectral densities and Gi(t,x,y) and G2(t,x,y) 
the corresponding kernels for the operators g(tHi) and g(tH 2 ), respectively. Then 

G 1 {t,x,y) = G 2 {t,x,y)+o(t°°), in V', as t -> + , (110) 

in £ (Z7 x [/); and when g admits an extension that belongs to /C this also holds pointwise 
in t. 

Let us consider some illustrations. 

Example. Let K(t,x,y) = (e(x,y;X), e~ xt ) be the heat kernel, corresponding to the 
operator K(t) = e~ tH , so that 

OK 

— = -HK, t>0, (111) 

and 

K(0 + ,x iy ) = 5(x-y). (112) 
In this case g(t) = x(^) e_ * admits extensions in /C. Thus the expansions 
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. ~ (-l) n H n 5(x -y)t n , . , 

K(t,x,y)~J2— f — ast^O + (113) 

n=0 n - 

in the space Z>'(.M x M), and 

fsT(t, x, y) = o(t°°) as t -> + , with x ^ y, (114) 

hold pointwise in £. ■ 

Example. Let U(t,x,y) = (e(x,y;X), e~ %xt ) be the Schrodinger kernel, corresponding 
to I7(t) = e~ itH , so that 

i— = HU, t>0 (115) 

and 

[/(0 + ,x,y)= ( 5(x-i/). (116) 
Here the function e~ lt belongs to S' but not to /C. Therefore, the expansions 

s ™(-i) n H n 5(x-y)t n , , . 

U(t,x,y)~J2— — , ast^0+, in I?', (117) 

n=0 n - 

and 

U(t,x,y) = o(0, in:p/ ' ast^0 + with x ^ y, (118) 

are distributional or "averaged" in t. 

Consider, for instance, U(t,x,y) = (A7rt)~ n ^ 2 e~^ x ~ y ^ Z 4 *, corresponding to the case when 
H is the negative Laplacian, —A, acting on M. = R n . If x ^ y are fixed, U (t, x, y) oscillates 
as t — > + , but ( |118| ) holds in the distributional sense. (This also follows from the results of 
Ref. m, because lim e ias = (C) if a ^ 0.) 

Notice that it also follows that if U% is the Schrodinger kernel associated to —A in any 
open subset M. of R n with boundary conditions that make the operator positive and self- 
adjoint, then 

Ux(t,x,y) = {A7Tt)- n / 2 e^ x - yl2/4t + o(t°°), in V, as t -> 0+ (119) 

uniformly and strongly on compacts of Ai. In particular, 

Ux{t,x,x) = (47rfr n/2 + o(t°°), inP', as t -> + , (120) 

uniformly and strongly on compacts of M.. (We reiterate that this holds in the interior of 
M only) ■ 

Apart from ( |120]) we have had little to say in this section about the much-studied diagonal 
expansions of the heat and Schrodinger kernels. See, however, Ref. |21, where the methods 
of this paper are applied in that arena. 
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VII. EXPANSION OF GREEN KERNELS II: GLOBAL EXPANSIONS 



When considering a second-order differential operator on a one- dimensional manifold, 
the variable A of the spectral density e(x, y\ A) is often replaced by the variable uj defined 
by uj 2 = A. For instance, the asymptotic behavior of ui n = A^ 2 , the square root of the 
nth eigenvalue, has a more convenient form than that for A n in the case of regular Sturm- 
Liouville equations. But, does this change of variable have any effect on the expansion of 
the associated Green kernels? 

Consider, for instance, the behavior of the kernel 

2 00 

e(x, y\ A) = — sin no; sinm/<5(A — n 2 ), (121) 

which we studied before. Let uj 2 = A, so that 

00 

e(x, y; uj) = e(x, y; uj 2 ) = — V] sinnx sinm/ 8{ui — n). (122) 

71 n=l 

The behavior of e(x, y\ uj) at infinity can be obtained by studying the parametric behavior, 
i.e., e(x, y\ auS) as a — > 00. Letting e = 1/a, we are led to consider the development of 

00 

— V^sinnx sin ny(j){en) = (e(x,y;uj), (j){en)) (123) 

71 n=l 

as e — > + . The even moments of e coincide with those of e of half the order; i.e., 

(e(x,y;uj),uj 2k )=6^\x- y ). (124) 

But we also have to consider odd-order moments, such as 

00 

J(x,y) = (e(x,y;uj),uj) = ^2n smnx sinm/. (125) 

n=l 

The operator corresponding to J is not the derivative d/dx, but rather (d/dx)Q, where Q is 
a Hilbert transform. Thus, J is not a local operator: in general supp J(4>) is not contained 
in supp 4>. Thus, the expansion 



e(x, y- „») ~ «t ~ «JM _ jfe + ^ - f V) - . . . ( 126) 

a a 2 1o A 



has a nonlocal character: The expansion of (e(x, y; cruj), f(y)) as a — > 00 may have nonzero 
contributions outside of supp/. 

The change of variable u; 2 = A, which seems so innocent, has introduced a new phe- 
nomenon into the expansion of the spectral density: the appearance of nonlocal terms. This 
will also apply to the expansion of the corresponding Green kernels. Consider the expansion 
of the cylinder kernel ( A9b| ) 

T ^,v) = <(x _ t y)2 + tr (127) 
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We have 



t . . t 5"(x-y)t 2 t 3 

5{x-y) + — v ' + — H as t -> 0" 



7r((x — t/) 2 + t 2 ) 7r(x — |/) 2 2 7r(x — y) 3 

(128) 

which is a nonlocal expansion because the odd-order terms are nonlocal. In particular, the 
expansion does not vanish when x ^ y. 

Why? Why do the results of the previous section fail? The operator corresponding to 
( |127| ) is e~ tHl/2 , where H = —d 2 /dx 2 on 1R. Thus, we arrive at the same question: How 
does the change A i— > A 1 / 2 = uj affect our results? However, the answer is now clearer: if 1 / 2 
is a nonlocal pseudo differential operator, and J(x, y) is simply its kernel. The nonlocality 
of the odd-order moments just reflects the fact that the basic operator H l l 2 is not local. 

The behavior of summability after changes of variables was already studied at the be- 
ginning of the century by HardytiE and is the central theme in Ref. [10| . 

When using the distributional approach, we can see that the change uj = A 1 / 2 , and similar 
ones, do not introduce problems at infinity. Rather, the point is that the change introduces 
a new structure at the origin. 

Let / 6 T>' have supp / C (0, oo). Then f{uj 2 ) is a well-defined distribution, given by 

(f(cj 2 ), 0H> = i</ (A), A" 1 / 2 0(A 1 / 2 )). (129) 

When G supp/, however, there is no canonical way to define f{uJ 2 )- That there are no 
problems at infinity follows from the results of Ref. |16| . 

Lemma 7.1. Let f be in V(M) with supp / C (0, oo). Then f is distributionally small 
at infinity if and only if fiyj 2 ) is. 

Proof: The generalized function / is distributionally small at infinity if and only if it 
belongs to K! . Thus it suffices to see that /(A) belongs to K! if and only if f{uj 2 ) does, and, 
by duality, it suffices to see that if supp0 C (0, oo) then </>(w) belongs to fC if and only if 
A- 1 / 2 0(A 1 / 2 ) does. ■ 

Therefore, if e(x, y;X) is the spectral density kernel corresponding to a positive self- 
adjoint operator H, then e(x,y;uj 2 ) is also distributionally small as a function of u, both 
distributionally and pointwise on x ^ y. However, as we have seen, the corresponding 
moment expansion for e(x,y;uj 2 ) will contain extra terms. These arise as a special case of 
a general theorem that extends the conclusions of Sections [V] and to the situation where 
the function or distribution g does not have a Taylor expansion at the origin. 

The spaces A{4> n } associated to an asymptotic sequence are discussed in Refs. |2^J2^. In 



particular, if a n is a sequence with 9fte a n /* oo, then the space IC{x an } consists of those 
smooth functions <f> defined on (0, oo) that have the behavior of the space K, at infinity but 
at the origin can be developed in a strong expansion 

<f>(x) ~ a lX ai + a 2 x a2 + a 3 x a3 + ■■■ as x -> 0+ (130) 

The point is that the a n need not be nonnegative integers. 
The functionals Sj € JC'{x a "} given by 

(6 j (x),(f>(x)) = a j (131) 
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play the role of the traditional delta functions. Each / G JC'{x a ™} admits a generalized 
moment asymptotic expansion, 

f(Xx) ~ f) K ^ X) asA-oo, (132) 

3=1 A 3 

where //(%) = (f(x),x aj ) are the moments. Therefore, if g G /C{x a ™} then the expansion 
of G(t) = (f(x),g(tx)) can be obtained from Ql32j ) as 

oo 

G(t) ~ J! ^ as * 0+ > ( 133 ) 
i=i 

where aj = (Sj(x), g(x)). But following the ideas of Section |V| we can define G(t) = 
(f(x),g(tx)) when g is a distribution of S' with suppg C [0, oo), whose behavior at the 
origin is of the form 

oo 

g{ex) ~ ^2 a j £aj x + as e — > + , (134a) 

3=1 

a fact that we express by saying that 

oo 

g(x) ~ 2J a i xQj 5 distributionally, as a; — >• + . (134b) 
i=i 

Then G(t) will have the same expansion (|133 ), but in the average or distributional sense. A 
corresponding result for operators also holds. 

Theorem 7.1. Let H be a positive self-adjoint operator on the domain X of the Hilbert 
space TC. Let be the intersection of the domains of H n for n G IN. Let g G S' with 
supp g C [0, oo) have a distributional expansion of the type g(x) ~ Y^jLi a j xC * j as x + ; 
where 3?ea n /* oo. Then G(t) = g(tH) can be defined as an element of S' \R, Li^X^^TL)) 
with support contained in [0, oo), and G(t) admits the distributional expansion 

oo 

G(t) ~ a 3 H ° 3 taj ' as t ~* 0+ ' m v '- ( 135 ) 

3=1 



When g belongs to JC{x an }, ( 135 ) becomes a pointwise expansion. ■ 

Therefore, we may generalize our previous discussion as follows: If H is a differential 
operator, then the expansion of the Green function of g(tH) is local or global depending on 
whether the expansion ( 134\) of g at the origin is of the Taylor- series type or not. 

Example. The small-t expansion of the cylinder function T(t,x,y) described in ( |A3| ) is 
given by 

r(«,,, g )~f (=£53ffe - y))t " a si ^0+. (136) 

n=0 U - 

The expansion is pointwise in t and distributional in (x,y). The expansion is also pointwise 
in Hor x y, but we do not get T(t, x, y) = o(t°°) because the odd terms in the expansion do 
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not vanish for x 7^ y. This type of behavior is typical of harmonic functions near boundaries. 
■ 

We may look at the locality problem from a different perspective. Suppose Hi and H2 
are two different self-adjoint extensions of the same differential operator on a subset U of 
M.. Then the two cylinder kernels T\{b,x,y) and T 2 (t,x,y) have different expansions as 
t — ► + even if (x, y) G U x U. The same is true of the associated Wightman functions. In 
both cases the small-t expansion of the Green kernel reflects some global properties of the 
operators H\ and H2 ■ Since the cylinder and Wightman functions are constructed from the 
operators Hi and H 2 , one may ask if this nonlocal character can already be observed in 
the spectral densities e H i/i(x,y\ A), j = 1,2. Interestingly, the nonlocal character cannot be 

seen in the Cesaro behavior, since according to Theorem 7.2 below we have 

e„i/2(x, y; A) = e„i/a(x, y; A) + o(A~°°) (C) as A — > 00 , (137) 

for (x, y) G C/ x U. Instead, the nonlocal character of the small-t expansion of the Green 
kernels is explained by the difference in the moments. (Recall Theorem 2.1 and the formula 



(0)0 

We finish by giving a result that justifies ( |137| ) and also has an interest of its own. 

Theorem 7.2. Let Hi and H 2 be two pseudodifferential operators acting on the manifold 
M., with spectral densities ej(x,y;X) for j = 1,2. Let U be an open set of M. and suppose 
that Hi — H 2 is a smoothing operator in U. Then 

ei(x, y; A) = e 2 (x, y; A) + o(A~°°) (C) as A -> 00, (138) 

in the topology of the space £{U x U) and, in particular, pointwise on (x,y) G U x U. 

Proof: If G T>(R), then 4>(Hi) — 4>(H 2 ) is a smoothing operator, thus (ei(x,y;X) — 
e 2 (x,y; A), f(x)g(y)) is a well-defined element of V(R) given by 

((ei(x, y; A) - e 2 (x, y; A), f{x)g(y)), 0(A)) = ((<j>(Hi) - </>(H 2 ))f, g). (139) 

In general this generalized function is not a distributionally small function of A, but if 
supp / <ZU and supp g C U then all the moments 

« ei (x, y; A) - e 2 (x, y; A), /(x)^(y)}, A") = ((fl? - i? 2 n )/, ^) (140) 

exist because Hi — H 2 is smoothing in {7. Therefore (ei(x, y; A) — e 2 (a;,y;A), f(x)g(y)) 
belongs to /C'(R); that is, it is a distributionally small function. Hence, 

(e 1 (a;,y;A)-e 2 (x,y;A) ) f(x)g(y)) = o(\-°°) (C) as A ^ 00 (141) 



for each /, g G £'(U), and (|138|) follows by duality 
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APPENDIX: THE SIMPLEST ONE-DIMENSIONAL EXAMPLES 



Let if be a positive, self-adjoint, second-order linear differential operator on scalar func- 
tions, on a manifold or region M.. We are concerned with distributions of the type G(t, x, y) 
(tell, x £ Ai, y £ Ai) that are integral kernels of parametrized operator- valued functions 
of H. In particular: 

(1) The heat kernel, K(t, x, y), represents the operator e~ tH , which solves the heat equa- 
tion 

-™=H% \ha*{t,x) = f{x), (Al) 

for (t, x) £ (0, oo) x M, by W(t, x) = e~ tH f(x). 

(2) The Schrddinger propagator, U(t,x,y), is the kernel of e~ %tH , which solves the 
Schrodinger equation 

i— = HS>, Km*(t,x) = f{x), (A2) 

for (t,x) £ R x M, by W(t,x) = e~ itH f(x). 

(3) Let T{t, x, y) be the integral kernel of the operator e~ 1 ^ , which solves the elliptic 
equation 

-—- + H^ = 0, lim*(t,a;) = f(x), lim tf(t, x) = 0, (A3) 

in the infinite half-cylinder (0, oo) x Ai, by ^(t,x) = e _iv/ ^/(x). We shall call this the 
cylinder kernel of H. It may also be regarded as the heat kernel of the first-order pseudo- 
differential operator \fH. 

(4) The Wightman function, W(t,x,y), is the kernel of (2\nT)~ 1 e~ lty/ ~^ '. This operator 
solves the wave equation 

d 2 ^ 

~w = B9 < A4 > 

with the nonlocal initial data 

\im^(t,x) = (2VH)- 1 f(x), \im—(t,x) = --f(x). (A5) 

The significance of W is that it is the two-point vacuum expectation value of a quantized 
scalar field satisfying the time-independent, linear field equation (|A4): 

W(t,x,y) = (0\<P(t,x)<P(0,y)\0) 

(e.g., Ref. ||, Chapters 3-5). 

These four kernels are rather diverse in their asymptotic behavior as t approaches 
and also in the convergence properties of their spectral expansions. Most of the relevant 
mathematical phenomena that distinguish them can be demonstrated already in the simplest 

case, 
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with M. either R or a bounded interval (0, tc). (For other one- dimensional Ai see Ref. |TU[ ) 
We record here the spectral expansion (Fourier transform or series) of each kernel and also 
its actual functional value (closed form or image sum). 



Heat kernel: 



Case M = R. 



K(t, x,y) = ^- [°° e ik{x - y) e- k2t dk. (A7a) 



2tt 



oo 



K(t, x, y) = (Ant)- ll2 e~ {x ~ y)2/u . (A7b) 



Schrddinger propagator: 



1 r 00 9 
U(t }X} y) = — e tk{x - y) e- lkt dk. (A8a) 

Z7T J-oo 

U(t, x, y) = e ~ l( - ssnt)n/4 (4irt)- l/2 e l{x ~ y)2/At . (A8b) 



Cylinder kernel: 



T(t, x,y) = ^- ^-v) e -\k\t dL ( A9a ) 

Z7T J-cxd 



Wightman function: 



T(t, x,y) = - ^— - . (A9b) 



1 Z" 00 ,e" i|fc| * 
W(t , x,y) = — e lk(x ~ v) —— dk. (A10) 

47T J-oo fc 



This integral is divergent at k = and does not make sense even as a distribution except on 
a restricted class of test functions. This "infrared" problem, which is irrelevant to the main 
issues of the present paper, disappears when one (1) goes to higher dimension, (2) adds a 
positive constant ("mass") to H, or (3) takes one or more derivatives of W with respect 
to any of its variables. Therefore, we do not list an integrated form of (|A10|) . (If we gave 



one, it would be nonunique and would grow logarithmically in (x — y) 2 , thus being useless 
for forming an image sum for (|A14j) .) For more information about infrared complications in 



simple model quantum field theories, see Ref. HO and references therein. 
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Case M = (0,vr) 



Heat kernel: 



K(t, x, y) = — ^2 sin(fcc) sin(A;y)e" 



-k 2 t 



7T 



k=l 



K(t,X,y) = (A7rty 1 / 2 £ ^-y-2Nn) 2 /4t _ e - {x+ y-2N^ /4t 

Schrddinger propagator: 



N=-oo 



U (t, x,y) = sin(A;a;) sin(%)e" 



-ikH 



TT 



k=l 



JJ^ ^ = e -j(sgnt)7r/4 (4^) -1/2 ^ ^ e i(x-y-2Nn) 2 /At _ e i(x+y-2N7r) 2 /At 

Cylinder kernel: 



N=-oo 



Z x -> 

T(t, x, y) = — B\n{kx) s\n{ky)e 



-kt 



TT 



k=l 



T(t,x,y) = - 



TT 



N=-oo 



(x-y- 2Nir) 2 + t 2 (x + y- 2Ntt) 2 + t 2 



In this case a closed form is obtainable: 



T(t,x,y) = ± 



sinht 



sinht 



cosh t — cos(x — y) cosh t — cos(x + y) 



Wightman function: 



W(t, x, y) = ^2 sin(fcc) sin(%) 



-ikt 



TT 



k=l 



k 



This is expressible in the closed form 



W(t,x,y) = ^hx 



cost — cos(a; + y) 



cost — cos(x — y) 



+ -P(t,x,y), 



where 



P(t,x,y) = < 



' —1 for 2kn — f(x + y) <t < 2krr — \x — y\, 

for 2kn — \x — y\ < t < 2kn + \x — y\, 

1 for 2kir + \x - y\ < t < 2kir + f(x + y), 

k for 2kir + f(x + y) <t < 2(k + 1)tt - f(x + y), 
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for fcG2i, with 



/(*) 



z if < Z < 71", 

27T — Z if 7T < Z < 27T. 



(A16) 



(P is essentially the standard Green function for the Dirichlet problem for the wave equation 
- the d'Alembert solution modified by reflections.) 

Proof of ( CT : Write flMSap as 



-| oo 



k(x+y)~kt ^ik(x—y)~kt ^—ik(x—y)~kt 



Evaluate each sum by the geometric series 



fe=i 



kz 



l-e z 



(A17) 



Some algebraic reduction yields ( |A13c 



Proof of ( \A14ty : Start with the well-known dispersion relation 



1 ^1 _ T, \ 

— — - = P~T mb(x) 
x ± x 



(A18) 



where Vj- is a principal- value distribution, whose antiderivative is In \x\ + constant. This 
generalizes easily to: 

Lemma A.l. Let F be analytic in a region Q \ {xi, . . . , x n }, Xj G M, where Q intersects 
M on (a, h) and each Xj is a simple pole of F with residue ctj . Then 



F(x - iO) = VF(x) + m aj5(x - 



x 



if 



(A19) 



To apply this to ( |A14b| ), replace t by ib in ( |A13a| ) and (|A13c|) : 

^ . /, x . /7 x -ikb i ( sin ^ sin ^ 
> sm fa sm(ky)e 1 — - ; 7 — 

fc=i 



4 V cos 6 — cos(x — y) cos 6 — cos(x + y) 



for Im 6 < 0. There are poles at b = 2kir ± (x — y) with residue if 4 and at 6 = 2/ot ± (x + y) 
with residue —i/4. Thus for 6 6 R, 



sin fa sin fcy e 



-ikb 



k=l 



V 



sin 6 



sin 6 



cos 6 — cos(x — y) cos 6 — cos(x + y) 



00 

+ ni Y - 2kn - x + y) + S(b - 2krc + x - y) - S(b - 2kn - x - y) - 5(b- 2kn + x + 

fc=— 00 

(A20) 
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Also, integrating ( |A13a|) -( |A13c|) and letting t — > + , one gets 



sin kx sin ky 1 / 1 — cos(x + y) 
fr[ k 4 \1 - cos(x - y) 



(A21) 



Therefore, integrating ( |A20|) yields 



E 

*;=i 



sin fcx sin ky e lkt 
k 



„l oo oo 

/ ^ sin kx sin % e~ ifcfe db + ^ 
fe=i fe=i 



sin fcx sin % 



4 jo ycoso — cos(x — y) coso — cos(x + y) J 4 \1 — cos(x — 



+ — / V (<5(6-2A;7r-x + ?/) + 5(6-2A;7r + x-y)-5(6-2A;7r-x-y)-5(6-2A;7r + x + ?/)) db 
4 jo ' 



fc=— oo 



cost — cos(a; + y) 



cost — cos(x — y) 



+ jP(t,x,y), 



which is ( |A14b| ). 
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